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Abstract 

First observations and measurements of the branching fractions of the B® — > D + D~ , 
B® — > D + D~ and B Q S — > D°D° decays are presented using 1.0 fb _1 of data collected 
by the LHCb experiment. These branching fractions are normalized to those of 
B° — > D + D~, B° — > D + D~ and B~ — > D°D~, respectively. An excess of events 
consistent with the decay B° — > D°D° is also seen, and its branching fraction is 
measured relative to that of B~ — > D°Dj . Improved measurements of the branching 
fractions B(B® — > DfD~) and B{B~ — > D°D~) are reported, each relative to 
B(B° — > D + D~). The ratios of branching fractions are 



1.08 ±0.20 ±0.10, 
0.050 ±0.008 ±0.004, 
0.019 ±0.003 ±0.003, 
0.0014 ±0.0006 ±0.0002, 
0.56 ±0.03 ±0.04, 
1.22 ±0.02 ±0.07, 
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where the uncertainties are statistical and systematic, respectively. 
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Double-charm decays of B mesons can be used to probe the Cabibbo-Kobayashi-Maskawa 
matrix [lj[2] elements, and provide a laboratory to study final state interactions. The 
time-dependent CP asymmetry in the B° — > D + D~ decay provides a way to measure the 
B° mixing phase |3j[4], where information from other double-charm final states can be 
used to account for loop (penguin) contributions and other non-factorizable effects [H]-[9]. 
Double-charm deca ys of B mesons can also be used to measure the weak phase 7, assuming 

The purely CP-even B° s 



[/-spin symmetry 
using the B® 



10 



11 



— > DfD~ decay is also of interest, 
as it can be used to measure the B° s mixing phase. Moreover, a lifetime measurement 
— > D+D~ decay provides complementary information on Ar s 11 13] to 
that obtained from direct measurements 
eigenstates 



15 16 



14 , or from lifetime measurements in other CP 



The study of B DD' decayeQ can also provide a better theoretical understanding 
of the processes that contribute to B meson decay. Feynman diagrams contributing to 
the decays considered in this paper are shown in Fig. Ml The B° s D°D°, B° s D + D~ 
and B° — > D°D° decays are mediated by the IV-exchange amplitude, along with penguin- 
annihilation contributions and rescattering 17 . The only other observed B meson decays 
of this type are B° £>i* )4 'K^~ and Bl 



of 10' 



18 



— > 7r + 7r , with branching fractions of the order 



and 10 -6 (l9j, respectively. Predictions of the B° s — » D + D~ branching fraction 
using perturbative approaches yield 3.6 x 10 -3 [20 , while the use of non-perturbative 
approaches has led to a smaller value of 1 x 10 -3 21 . More recent phenomenological 



studies, which assume a dominant contribution from rescattering, predict a significantly 
lower branching fraction of B(B° S D + D~) = B(B° S D°D°) = (7.8 ± 4.7) x 10" 5 [Tr] . 

This paper reports the first observations of the B° s — > D + D~, B° s — > D+D~ and 
B° s — > D°D° decays, and measurements of their branching fractions normalized relative to 
those of -B — > D + D~ , B° — > D+D~ and B~ — > D°D~ , respectively. An excess of events 
consistent with B° — > D°D° is also seen, and its branching fraction is reported. Improved 
measurements of the ratios of branching fractions B(B° S — > D+D~)/B(B° — > DfD~) and 
B(B~ — > D°D~)/B(B° — > D+D~) are also presented. All results are based upon an 
integrated luminosity of 1.0 fb^ 1 of pp collision data at y/s = 7 TeV recorded by the LHCb 
experiment in 2011. Inclusion of charge conjugate final states is implied throughout. 



32 2 Data sample and candidate selection 

33 The LHCb detector [22 1 is a single-arm forward spectrometer covering the pseudorapidity 

34 range 2 < r] < 5, designed for the study of particles containing b or c quarks. The detector 

35 includes a high precision tracking system consisting of a silicon-strip vertex detector 

36 surrounding the pp interaction region, a large-area silicon-strip detector located upstream 

37 of a dipole magnet with a bending power of about 4 Tm, and three stations of silicon-strip 

throughout this paper, the notation D is used to refer to a D + , D° or D+ meson, and B represents 
either a B°, B~ or 5° meson. 



1 




Figure 1: Feynman diagrams contributing to the double-charm final states discussed in this 
paper. They include (a) tree, (b) TF-exchange and (c) penguin diagrams. 



38 detectors and straw drift tubes placed downstream. The combined tracking system has a 

39 momentum resolution (Ap/p) that varies from 0.4% at 5GeV/c to 0.6% at lOOGeV/c, and 
an impact parameter (IP) resolution of 20 |xm for tracks with high transverse momentum 
(Pt)- The impact parameter is defined as the distance of closest approach of a given 
particle to the primary pp interaction vertex (PV). Charged hadrons are identified using 



two ring-imaging Cherenkov detectors 23 . Photons, electrons and charged particles are 



44 identified by a calorimeter system consisting of scintillating-pad and preshower detectors, 

45 an electromagnetic calorimeter and a hadronic calorimeter. Muons are identified by a 

46 system composed of alternating layers of iron and multiwire proportional chambers. 



The trigger 24 consists of a hardware stage, based on information from the calorimeter 
and muon systems, followed by a software stage that performs a partial event reconstruction 
(only tracks with > 0.5 GeV/c are reconstructed and used). The software trigger 
requires a two-, three- or four-track secondary vertex with a large track j?t sum and a 

51 significant displacement from any of the reconstructed PVs. At least one track must have 

52 > 1.7GeV/c and IP x 2 greater than 16 with respect to all PVs. The IP x 2 is defined 

53 as the difference between the x 2 °f the PV reconstructed with and without the considered 



54 particle. A multivariate algorithm 25 is used to identify secondary vertices that originate 

55 from the decays of b hadrons. 

56 For the ratios of branching fractions between modes with identical final states, no 

57 requirements are made on the hardware trigger decision. When the final states differ, a 

58 trigger selection is applied to facilitate the determination of the relative trigger efficiency. 

59 The selection requires that either (i) at least one of the tracks from the reconstructed 
eo signal decay is associated with energy depositions in the calorimeters that passed the 
6i hardware trigger requirements, or (ii) the event triggered independently of the signal decay 
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62 particles, e.g., on the decay products of the other b hadron in the event. Events that do 



63 



103 



not fall into either of these two categories (~5%) are discarded. 



64 Signal efficiencies and specific backgrounds are studied using simulated events. Proton- 



65 proton collisions are generated using Pythia 6.4 26 with a specific LHCb configura- 
ee tion 27 . Decays of hadronic particles are described by EvtGen 28 in which final state 
67 radiation is generated using Photos |29|. The interaction of the generated particles 

as 



es with the detector and its response are implemented using the Geant4 toolkit 30 

eg described in Ref. [31 1 . Efficiencies for identifying K + and ?t + mesons are determined using 

70 D* + calibration data, with kinematic quantities reweighted to match those of the signal 

71 particles |23j . 

72 Signal B candidates are formed by combining pairs of D meson candidates recon- 

73 structed in the following decay modes: D° — > K~n + or K~7t + 7t~tt + , D + — > K~7t + tt + and 

74 D+ K+K~ir + . The D° -> K-tt+it-tt + decay is only used for B° (s) D°D° candidates, 

75 where a single D° — > K~7t + 7t~7t + decay in the final state is allowed, which approximately 

76 doubles the total signal efficiency. A refit of signal candidates with D mass and vertex 

77 constraints is performed to improve the B mass resolution. 

78 Due to similar kinematics of the D + — > K~h + 'k + . Df — > K + K~ix + and — > pK~7i + 

79 decays, there is cross-feed between various 6-hadron decays that have two charm particles 
in the final state. Cross- feed between D + and occurs when the K~n + h + invariant 
mass is within 25 MeV/c 2 (~ 3 times the experimental resolution) of both the D + and 
-D+ masses under the h + = n + and h + = K + hypotheses, respectively. In such cases, 
an arbitration is performed as follows: if either \M(K + K~) — < 10 MeV/c 2 or h + 

84 satisfies a stringent kaon particle identification (PID) requirement, the D candidate is 

85 assigned to be a Df meson. Conversely, if h + passes a stringent pion PID requirement, 
se the D candidate is taken to be a D + meson. Candidates that do not pass either of these 

87 selections are rejected. A similar veto is applied to D + and decays that are consistent 

88 with the — > pK~Tr + decay hypothesis if the proton is misidentified as a ?t + or K + , 

89 respectively. The efficiencies of these D selections are determined using simulated signal 

90 decays to model the kinematics of the decay and D* + — > D°tt + calibration data for the 

91 PID efficiencies. Their values are given in Table [TJ 

92 To suppress contributions from non-DD' final states, the reconstructed D decay vertex 

93 is required to be downstream of the reconstructed B decay vertex, and the B and D decay 

94 vertices are required to have a vertex separation (VS) x 2 larger than two. Here, the VS y 2 

95 is the difference in y 2 between the nominal vertex fit and a vertex fit where the D is 

96 assumed to have zero lifetime. The efficiencies of this set of requirements are obtained 

97 from simulation and are included in Table [TJ 

98 To further improve the purity of the B — > DD' samples, a boosted decision tree 



99 (BDT) discriminant is used to distinguish signal D mesons from backgrounds 32 33 . The 

100 BDT uses five variables for the D meson and 23 for each of its children. The variables 

101 include kinematic quantities, track quality, and vertex and PID information. The signal 

102 and background distributions used to train the BDT are obtained from B° — > D + n~ , 
B~ — > D°7T~ and B° s — > Dfir~ decays from data. The signal distributions are background 



104 subtracted using weights [34] obtained from a fit to the B candidate invariant mass 
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Table 1: Individual contributions to the efficiency for selecting the various B — > DD' final states. 
Shown are the efficiencies to reconstruct and trigger on the final state, and to pass the charm 
cross- feed veto, the VS % 2 and BDT selection requirements. The total selection efficiency is the 
product of these four values. The relative uncertainty on the selection efficiency for each decay 
mode due to the finite simulation samples sizes is 2%. Entries with a dash indicate that the 
efficiency factor is not applicable. 



B° s -+ D+D- 



B° -»■ D+D (loose selection) 



B\ s) -»• D°D°, 



7i 



,K + 7T- 



B~ -> D°D7 



Rec.xTrig. 

0.140 
0.130 
0.447 
0.128 
0.238 



Efficiencies (%) 
Cross-feed veto 

88.4 
77.8 



92.5 



VS X 2 

75.4 
82.9 
73.7 
74.6 
75.0 



BDT 

97.5 
100.0 
57.8 
63.6 
99.2 



105 distribution. The background distributions are taken from the high B mass sidebands in 

we the same data sample. 

107 It is found that making a requirement on the product of the two D meson BDT responses 

loa provides better discrimination than applying one to each BDT response individually. The 

109 optimal BDT requirement in each decay is chosen by maximizing N$ / a/ N$ + N-q. The 

no number of signal events, iVg, is computed using the known (or estimated, if unknown) 

in branching fractions, selection efficiencies from simulated events, and the BDT efficiencies 

112 from the B° — > D + tt~ , B~ — > D°tt~ and B® — > Dfn~ calibration samples, reweighted to 

in account for small differences in kinematics between the calibration and signal samples, 

lu The number, is the expected background yield for a given BDT requirement. The 

115 efficiencies associated with the optimal BDT cut values, determined from an independent 

lie subset of the B — > Dn~ data, are listed in Table [I] Correlations between the BDT values 

117 for the two D mesons are taken into account. 

us For the purpose of measuring B(B Q S — > D+D~)/B(B° — >■ D+D~), only loose BDT 

n9 requirements are imposed since the expected yields are relatively large. On the other hand, 

120 for B(B° S -»■ D+D-)/B(B° D+D~), the expected signal yield of B° s -»■ D+D~ decays is 

121 small; in this case both the signal and normalization modes are required to pass the same 

122 tighter BDT requirement. The different BDT selections applied to the B° — > D+D~ decay 

123 are referred to as the "loose selection" and the "tight selection." Since the final state 

124 is identical for the tight selection, the BDT efficiency cancels in the ratio of branching 

125 fractions, and is not included in Table [TJ 

126 For B ,, — > D°D° candidates, a peaking background from B — > D* + n~ — > (D°7t + )tt~ 

127 decays, where the ir + is misidentified as a K + , is observed. This contribution is removed 

128 by requiring the mass difference, M(K~tt + tt + ) — M(K~tt + ) > 150 MeV/c 2 , where the K + 

129 in the reconstructed decay is taken to be a 7r + . After the final selection around 2% of 
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130 events in the B° s — > DfD s decay mode contain multiple candidates; for all other modes 

131 the multiple candidate rate is below 1%. All candidates are kept for the final analysis. 



3 Signal and background shapes 

The B — > DD' signal shapes are all similar after the D mass and vertex constraints. The 



signal shape is parameterized as the sum of two Crystal Ball (CB) functions 35 , which 
account for non-Gaussian tails on both sides of the signal peak. The asymmetric shapes 
account for both non-Gaussian mass resolution effects (on both sides) and energy loss 
due to final state radiation. The two CB shapes are constrained to have equal area and 
a common mean. Separate sets of shape parameters are determined for B° — > D+D~, 
B° s — > DfD~ and B~ — > D°Dj using simulated signal decays. In the fits to data, the 
signal shape parameters are fixed to the simulated values, except for a smearing factor 
that is added in quadrature to the widths from simulation. This number is allowed to 
vary independently in each fit, but is consistent with about 4.6 MeV/c 2 across all modes, 
resulting in a mass resolution of about 9 MeV/c 2 . For the more rare £?9> — > D°D° and 

B® s ^ — > D + D~ decay modes, the B° s — > D+D~ signal shape parameters are used. In 

determining the signal significances, the signal shape is fixed to that for B° s — > D+D~, 
including an additional smearing of 4.6 MeV/c 2 . The impact of using the B° — > D+D~ or 
B~ — > D°D~ signal shapes on the signal significances is negligible. 

Several specific backgrounds contribute to the DD' mass spectra. In particular, decays 
such as B — > D^D*, where the D* mesons decay through pion or photon emission, produce 
distinct structures in all decays under consideration. The shapes of these backgrounds 
are derived from simulation, which are corrected for known resolution differences between 
data and simulated events, and then fixed in fits to the data. The relative yield of the two 
peaks in the characteristic structure from the decay D* — > D°tt is allowed to vary freely, to 
enable better modeling of the background in the low mass region. Since this mass region 
is significantly below the signal peaks, the impact on the signal yield determinations is 
negligible. 

A source of peaking background that contributes to B — > DDf modes are the 
B DK*°K + DK~ir + K + decays, where the K*°K + is not produced in a Df de- 
cay. Although the branching fractions for these decays [36] are about twice as large as 
that of the B — > DDf — > DK + K~n + decay channel, the 25 MeV/c 2 mass window around 
the known Df mass and the VS x 2 > 2 requirement reduce this contribution to about 1% 
of the signal yield. This expectation is corroborated by studying the Df candidate mass 
sidebands. The shape of this background is obtained from simulation, and is described 
by a single Gaussian function which has a width about 2.5 times larger than that of the 
signal decay and peaks at the nominal B meson mass. 

After the charm cross-feed vetoes (see Sect. [2]), the cross-feed rate from B° — > DfD~ 
decays into the B° s — > DfD~ sample is (0.7 ± 0.2)%. The shape of this misidentification 
background is obtained from simulation. A similar cross-feed background contribution 
from A° — > A+D~ decays is also expected due to events passing the A£ veto. Taking 
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Figure 2: Invariant mass distributions for (left) B® — > D+Dj and (right) B° — >• DfD~ 
candidates in the data with the loose BDT selection applied to the latter. The signal and 
background components are indicated in the legend. The A® — > A+Dj, B® — > Df K~ K + n~ and 
B° — > D~ K + K~tt + background components are too small to be seen, and are excluded from 
the legends. 



170 into account the observed yields of these decays in data, we fix the B — > D^D~ and 

171 A° — > A+D~ cross-feed yields to 35 and 15 events, respectively. Investigation of the D 

172 mass sidebands reveals no additional contributions from non-DD' backgrounds. 

173 The combinatorial background shape is described by an exponential function whose 

174 slope is determined from wrong-sign candidates. Wrong-sign candidates include the D+D+, 

175 D°D°, or D°(K + n~)D~ final states, in which no signal excesses should be present (ne- 

176 glecting the small contribution from the doubly Cabibbo suppressed B~ — > D°(K + tt~)D~ 

177 decay). For the B9 s — > D + D~ decay, the exponential shape parameter is allowed to vary 

178 in the fit due to an insufficient number of wrong-sign D + D + candidates. 



179 4 Fit results 

wo Figure [2] shows the invariant mass spectra for B° s — > DfD~ and B° — > D+D~ candidates. 

181 The results of unbinned extended maximum likelihood fits to the distributions are overlaid 

182 with the signal and background components indicated in the legends. Signal yields of 

183 451 ± 23 £° ->• D+D- and 5157 ± 64 B° ->• D+D~ decays are observed. 

184 Figure shows the invariant mass spectrum for B° — > DfD~ and B® — > DfD~ 

185 candidates, where the tight BDT selection requirements have been applied as discussed 

186 previously. We observe 36 ± 6 B° s — > D+D~ signal decays, with 2832 ± 53 events in 

187 the B° — > DfD~ normalization mode. The statistical significance of the B° s — > DfD~ 

188 signal corresponds to 10a by computing a/— 2 \n(C / C max ) , where £ max and C are the fit 
las likelihoods with the signal yields allowed to vary and fixed to zero, respectively. Variations 

190 in the signal and background model have only a marginal impact on the signal significance. 

191 The B° s — > D~Df decay is thus observed for the first time. 

192 The invariant mass spectrum for 5? — > D + D~ candidates is shown in Fig. 
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Figure 3: Invariant mass distribution for B° — > DfD~ and B® — > DfD~ candidates in the 
data, with the tight BDT selection applied. The distribution is plotted on a (left) linear and 
(right) logarithmic scale to highlight the suppressed B® — > DfD~ signal. Signal and background 
components are indicated in the legend. 



Peaks are seen at both the B° and B® meson masses, with yields of 165± 13 and 43±7 signal 
events, respectively. In the lower mass region, two prominent peaks from B° — > D* + D~ and 
B° — > D + D*~ decays are also evident. The significance of the B° s — > D + D~ signal yield is 
computed as described above, and corresponds to 11 a, establishing the first observation of 
this decay mode. 

Figure [4] (right) shows the D°D° invariant mass distribution and the results of the fit. 
Both (K~tc + , K + n~) and (K~tt + , K + n~n + n~) combinations are included. AB°-> D°D° 
signal is seen with a significance of lla, which establishes the first observation of this 
decay mode. The data also show an excess of events at the B° mass. The significance of 
that excess corresponds to 2Aa, including both the statistical and systematic uncertainty. 
The fitted yields in the B° s -> D°D° and B° D°D° decay modes are 45 ± 8 and 
13 ± 6 events, respectively. If both the B° s — > D°D° and B° — > D°D° decays proceed 
through W -exchange diagrams, one would expect the signal yield in B° — > D°D° to be 
~ (fd/fs) x \V cd /V cs \ 2 ~ 0.2 of the yield in B° s ->• D°D°, where we have used \V cd /V cs \ 2 = 



0.054 [18] and f s /f d = 0.256 ± 0.020 37 . The fitted yields are consistent with this 
expectation. The decay B~ — > D°D~ is used as the normalization channel for both 
the B® — > D°D° and B° — > D°D° branching fraction measurements, where only the 
D° — > K~n + decay mode is used. The fitted invariant mass distribution for B~ — > D°Dj 
candidates is shown in Fig. |5j The fitted signal yield is 5152 ± 73 events. 

The measured yields, N B ^ D j=),, relevant for the branching fraction measurements are 
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Figure 4: Invariant mass distributions for (left) B9 — > D+D and (right) B9s — > D°D° 
candidates in the data. Signal and background components are indicated in the legend. 
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Figure 5: Invariant mass distribution for B~ — > D°D~ candidates in the data. Signal and 
background components are indicated in the legend. The B~ — > D°K~ K+ir~ background 
components are too small to be seen, and are excluded from the legend. 



213 summarized in Table [2] The branching fractions are related to the measured yields by 



B(B° a D+D[ 



B(B° ->■ D+D- 
B(B° S -». D+D 



B(B° D+D' 
B{B° S -)• D+D 



B{B° ->■ D+D^ 



->• D°D~ 



fd 
fs 

fd 
fs 

fd 
fs 
fd 
fs 

e rel 



• e 



B°/B° 



rel 



• Crel 



K 



N B°^D+D- 



N 



B°^D+D- 



N 



e re i • k ■ 



B°^D+D~ 



N 



B°^D+D~ 



N 



■ e 



B°^D°D° 



rel 



N 



B~ 



N 



N 



B~ 



N 



B{D+ -> K-ix+ix^ 
B(D+ K+K'7i+) N j 



bo^d+d: 



B°^DJD- 



(1) 

(2) 
(3) 
(4) 
(5) 



S 



Table 2: Summary of the observed signal and normalization mode yields and their relative 
efficiencies, as used in the measurements of the ratios of branching fractions. The quoted 
uncertainties are statistical only. 



Measurement 


Signal 


Norm. 


Rel. eff. 




yield 


yield 


fc rel 


B(B° -KD+ d-) 
B(B°^D+D-) 


451 ± 23 


5157 ±64 


0.928 ±0.027 


B(B a a ^D+D-) 
B{B°^DfD-) 


36 ±6 


2832 ± 53 


1.0 


B{B°-+D+D-) 
B(B ^D+D~) 


43 ±7 


165 ± 13 


1.0 


B(B° S ^D°D ) 
B{B~^D°D7) 


45 ±8 


5152 ± 73 


0.523 ±0.016 


B(B°^D°D ) 
B(B~^D°D7) 


13 ±6 


5152 ± 73 


0.523 ±0.016 


B(B~^D D~) 
B(B°^D+D-) 


5152 ± 73 


5157 ±64 


0.508 ±0.011 



B{B~ -» D°DJ) _ B o /B - B(D+ -» K-n+n+) N B ^ D0DJ 

B(B° D+D-) ' B{D° K-n+) ' N B0 ^ DJD . " 1 } 

2u Here, it is assumed that B~ and B° mesons are produced in equal numbers. The relative 

215 efficiencies, e re i, are given in Table [2] They account for geometric acceptance, detection and 

216 trigger efficiencies, and the additional VS x 2 i BDT, and charm cross-feed veto requirements. 

217 The first four of these relative efficiencies are obtained from simulation, and the last two 

218 are data-driven. The indicated uncertainties on the relative efficiencies are due only to the 

219 finite sizes of the simulated signal decays. The average selection efficiency for B~ — > D°Dj 

220 relative to B° (s) D°D° is 

, e B ^ D0D7 B(Df -> K+K-n+)B(D° -» K~n+) 

Crcl e KVtKv [B(D° -> K-k+)Y + 2e K ^ K „B(D» -> K~>k+)B{D* K-n+n-n+) ' 1 j 

221 where the quantities e B -^ D o D7 = (0.166 ± 0.003)%, e K „ }K * = (0.190 ± 0.003)% 

222 and eKnmT,Kir = (0.061 ± 0.002)% are the selection efficiencies for the B~ — > D°D~ , 

223 B° s -> (D '^ K-n + , D° -> K+-R-) and B° s -> (D° -> R-7V+, D° -> A^+tt-tt+tt-) de- 

224 cays, respectively. The .D branching fractions, <B(D° — > K~n + ) = (3.88 ± 0.05)%, 

225 B(D° ->■ A-7r+7r-7r+) = (8.07 ± 0.20)%, £(£>+ ->. K+R-TT+) = (5.49 ± 0.27)%, and 

226 £(£>+ -»■ A-7T+7T+) = (9.13 ± 0.19)% are taken from Ref. [18] . 

227 The factor k is a correction that accounts for the lower selection efficiency associated 

228 with the shorter-lifetime CP-even eigenstates of the B° s system compared to flavor-specific 

229 final states [l4]. The impact on the B° s acceptance is estimated by convolving an exponential 

230 distribution that has a 10% smaller lifetime than that in flavor-specific decays with the 
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231 simulated lifetime acceptance. The resulting correction is k = 1.058 ± 0.029. In the B c 



232 sector, AT d /T d is below 1% 38 , and the lifetime acceptance is well described by the 

233 simulation. 

234 The measured ratios of branching fractions are computed to be 

B{B° S ->■ D+D 



B(B° D+D 
B(B° S D+D' 



B(B° ->■ D+D- 



= 1.08 ± 0.20 (stat) ± 0.10 (syst), 

= 0.050 ± 0.008 (stat) ± 0.004 (syst), 

= 0.019 ± 0.003 (stat) ± 0.003 (syst), 

: 0.0014 ± 0.0006 (stat) ± 0.0002 (syst) 

< 0.0024 at 90% CL ], 

: 0.56 ± 0.03 (stat) ± 0.04 (syst), 

= 1.22 ± 0.02 (stat) ± 0.07 (syst). 



the results obtained using 
)D°(K+n-) final states differ by 



the 
less 



23 5 For B{B^ D°D°)/B{B- D°D~), 

236 D°(K-n+)D°(K+n-n+n-) and D°(K~ir~ 

237 than one standard deviation. For the B° — > D°D° decay, we provide both the central value 

238 and the 90% confidence level (CL) upper limit. The upper limit is obtained by convolving 

239 the fitted likelihood with a Gaussian function whose width is the total systematic error, 

240 and integrating over the physical region. 



241 5 Systematic uncertainties 

242 A number of systematic uncertainties contribute to the measurements of the ratios of 

243 branching fractions. The sources and their values are summarized in Table [3| The dominant 

244 source of uncertainty on the branching fraction ratios comes from the b fragmentation 

245 fraction ratio, fd/f s , which has a total uncertainty of 7.8% (37), of which 5.3% is from the 

246 ratio of branching fractions B(D+ — > K+ K~ tt+) j B(D+ — > K~7t+tt+). For clarity, we have 

247 removed that portion of the uncertainty from fd/f s , and included its contribution in the 

248 row labeled B(D) in Table § For B(B° S -»■ D+Dj)/B(B° -»■ D+D~), the above D+/D+ 

249 branching fraction ratio from fd/f s cancels with the corresponding inverted ratio in Eq. [TJ 

250 On the other hand, in the ratio B(B° (s) -»■ D°D°)/B(B- D°Dj), the D+ -»■ K+R-n+ 

251 branching fraction enters as the square, after considering the D branching fractions used 

252 in computing fd/fs (see Eq. [4]). As a result, the uncertainty from B(D+ — > K+K~ir+) 

253 contributes 9.8% to the total uncertainty on B(B° (s) -> D°D°)/B(B- ->■ D°Dj); smaller 
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contributions from the limited knowledge of B(D° — > K~n + ) [1.3%], B(D° — > K~7t + 7t~tt + ) 
[2.5%] and B(D + -)• K~tt + -k + ) [2.1%] are also included in the 13(D) uncertainties. 

Another significant uncertainty results from the precision on 6-hadron lifetimes and 
decays of B° and B® to CP eigenstates. Using the measured value of the width difference, 



Ar s = 0.116 ± 0.018 ± 0.006 ps 39 we conservatively assume the CP-even lifetime to 
be in the range from 0.85 to 0.95 times the flavor-specific decay lifetime. With this allowed 
range a 2.9% uncertainty on the efficiencies for B° s decays to CP eigenstates is found. The 
average B® lifetime is known only to a precision of 3%, which leads to a 1.5% uncertainty 
on the selection efficiencies for B° s decays to flavor-specific final states. The B° and B~ 
lifetimes are known with sufficient precision that the associated uncertainty is negligible. 

Several of the efficiency factors are estimated from simulation. Most, but not all, of 
the associated systematic uncertainties cancel due to the similar or identical final states 
for the signal and normalization modes. For modes with an unequal number of tracks 
in the final state, a 1% uncertainty due to small differences in the IP resolution between 
data and simulation is assigned. The efficiency of the VS x 2 requirement is checked 
using the large B° — > D+D~ signal in data, and the agreement to within 1% with the 
efficiency from simulation is the assigned uncertainty. For B(B~ — > D°D~)/B(B° — > 
D+D~), a 1% uncertainty is attributed to the efficiency of track reconstruction. For 
B(B° S D°D°)/B(B- D°Dj), the one fewer track in the D (K<k)D°(K-k) final state is 
offset by the one extra track in D c \Kit)D q '(Kttitit) , relative to D° (Kir)D~(KK7r), leading 
to a negligible tracking uncertainty. The mass resolution in data is slightly larger than 
in simulation, resulting in slightly different efficiencies for the reconstructed D°, D + and 
invariant masses to lie within 25 MeV/c 2 of their known masses. This introduces 
a maximum of 1% uncertainty on the relative branching fractions. To estimate the 
uncertainty on the trigger efficiencies determined from simulation, the hadron trigger 
efficiency ratios were also determined using data. These efficiencies were measured using 
trigger-unbiased samples of kaons and pions identified in D* + — > D°tt + decays. Using 
this alternative procedure, we find that the simulated trigger efficiency ratios have an 
uncertainty of 2%. The combined systematic uncertainties in the efficiencies obtained from 
simulation are given in Table [3] 

The limited sizes of the B — > Dtt~ calibration samples lead to uncertainties in the 
BDT efficiencies. The uncertainties on the ratios vary from 1.0% to 2.0%. The uncertainty 
on the efficiency of the _D( S ) and A£ vetoes is dominated by the PID efficiencies, but they 
only apply to the subset of D candidates that fall within the mass window of two charm 
hadrons, e.g., both the D + and -D+ mesons, which occurs about 20% of the time for Df 
decays. Taking this fraction and the uncertainty in the PID efficiency into account, the 
veto efficiencies are estimated to have uncertainties of 1.0% for the D + veto, 0.5% for the 
D+ veto, and 0.3% for the veto. 

The fit model is validated using simulated experiments, and is found to be unbiased. 
To assess the uncertainty due to the imperfect knowledge of the various parameters 
used in the fit model, a number of variations are investigated. The only non-negligible 
uncertainties are due to the B — > DK~ K + n~ background contribution, which is varied 
from 0% to 2%, and the cross-feed from B® — > DfD~ decays into the B° s — > D+D~ sample. 
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Table 3: Sources of systematic uncertainty and their values (in %) for the ratios of branching 
fractions of the indicated decays. For B(B° {s) -> D°D°)/B(B- -»■ D°Dj), the error on fd / f s 

only applies to the B® — > D°D° decay, as indicated by the values in parentheses. 



Source 




B° S ^D+D- 


B° S ^D+D~ 


B9 ,-s>D°D° 

W 


B-^D°DJ 


B°— > D~ 


B°—>Df D~ 


B°—>D+ D~ 


B~— lD° D~ 


B°—¥ Df D~ 


fd/fs 


5.7 


5.7 


5.7 


~ (5-7) 




B(D) 




5.3 


5.3 


10.2 


2.5 


B meson lifetimes 


2.9 


1.5 


2.9 


2.9 




Eff. from simulation 


2.4 






2.2 


2.6 


BDT selection 


1.4 






2.2 


1.4 


Cross-feed vetoes 


0.6 






0.5 


1.0 


D mass resolution 


1.0 






1.0 


1.0 


Fit model 


2.1 


0.5 


0.5 


1.7 


2.1 


Simulated sample size 


3.0 


3.0 


3.0 


3.0 


3.0 


Total 


8.0 


8.5 


8.9 


11.7(13.0) 


5.5 



297 The uncertainty varies from 1.7% to 2.1%. For B(B° S D + D-)/B(B° ->• D+D~) and 

298 B(B° S — > D+D~)/B(B° — > D+D~), we assign an uncertainty of 0.5%, which accounts for 

299 potentially small differences in the signal shape for B° and B° s decays (due to the B°-B® 

300 mass difference). Lastly, the finite size of the samples of simulated decays contributes 

301 3% uncertainty to all the measurements. In total, the systematic uncertainties on the 

302 branching fraction ratios range from 5.5% to 13.0%, as indicated in Table [3] 



303 6 Discussion and summary 

304 First observations and measurements of the relative branching fractions for the decays 

305 B° s — > D + D~ , B° s — > D+D~ and B® — > D°D° have been presented, along with measure- 

306 ments of B(B° S — > D+D~) and B{B~ — > D Q D~). Taking the world average values for 

307 B(B° D+D-) = (7.2 ± 0.8) x 10~ 3 [18], the absolute branching fractions are 



B[B- D°D;) = (8.6 ± 0.2 (stat) ± 0.4 (syst) ± 1.0 (norm)) x 10~ 3 , 
B(B° S D+DJ) = (4.0 ± 0.2 (stat) ± 0.3 (syst) ± 0.4 (norm)) x 10" 3 . 

The third uncertainty reflects the precision of the branching fraction for the normaliza- 
tion mode. These measurements are consistent with, and more precise than, both the 



current world average measurements [18 as well as the more recent measurement of 
B(B° S DfD~) (40]. 

The measured value of B(B° S D+D;)/B(B° D+D~) = 0.55 ± 0.06 is signif- 
icantly lower than the naive expectation of unity for the case that both decays are 
dominated by tree amplitudes (see Fig. [l](a)), assuming small non-factorizable effects and 
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comparable magnitudes of the B^ — > form factors |41|. Unlike B° — > D+D , the 

B° s — > D+D~ decay receives a contribution from the PU-exchange process (see Fig. [ijb)) , 

317 suggesting that this amplitude may not be negligible. Interestingly, when comparing the 

318 B° s — > D+D~ and B° — > D + D~ decays, which have the same set of amplitudes, one finds 

319 \V cd /V cs \ 2 -B(B° s ^DfD;)/B(B°^D+D-)~l. 
Using B(B° -> D+D~) = (2.11 ± 0.31) x 10~ 4 and B{B~ D°D~) =(10.0 ± 1.7) x 



316 



320 



321 10 3 [18], the following values for the branching fractions are obtained 

B(B° S D + D~) = (2.2 ± 0.4 (stat) ± 0.2 (syst) ± 0.3 (norm)) x 10" 4 , 
B(B° S ->• D°D°) = (1.9 ± 0.3 (stat) ± 0.3 (syst) ± 0.3 (norm)) x 10" 4 , 
B(B° D°D°) = (1.4 ± 0.6 (stat) ± 0.2 (syst) ± 0.2 (norm)) x 10" 5 . 



324 



328 



330 



20 



21 



322 The first of these results disfavors the predicted values for B(B® — > D + D~) in Refs. 

323 which are about 5-15 times larger than our measured value. The measured branching 
fractions are about a factor of 2-3 larger than the predictions obtained by assuming that 



325 these decay amplitudes are dominated by rescattering 17 . As discussed above for the 

326 B(B Q S — > D+D~) measurement, this may also suggest that the TU-exchange amplitude 

327 contribution is not negligible in B — > DD' decays. For precise quantitative comparisons of 
these B® branching fraction measurements to theoretical predictions, one should account 



329 for the different total widths of the CP-even and CP-odd final states 12 . 

The Cabibbo suppressed P° — > DfD~ decay is also observed for the first time. Its 



331 absolute branching fraction is 

B(B° S D+D~) = (3.6 ± 0.6 (stat) ± 0.3 (syst) ± 0.4 (norm)) x 10' 



332 This value is consistent with the expected suppression of | V c d/V cs \ 2 . 

333 The results reported here are based on an integrated luminosity of 1.0 fb -1 . A data 

334 sample with approximately 2.5 times larger yields in these modes has already been collected 

335 in 2012, and larger samples are anticipated in the next few years. These samples give 

336 good prospects for CP-violation measurements, lifetime studies, and obtaining a deeper 

337 understanding of the decay mechanisms that contribute to 6-hadron decays. 
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